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Multi-scaling properties of one-dimensional truncated Levy flights are studied. Due to the broken 
self-similarity of the distribution of jumps, they are expected to possess multi-scaling properties in 
contrast to the ordinary Levy flights. We argue this fact based on a smoothly truncated Levy 
distribution, and derive the functional form of the scaling exponents. Specifically, they exhibit 
bi-fractal behavior, which is the simplest case of multi-scaling. 



I. INTRODUCTION 

Levy flights are random processes based on Levy stable distributions [Q. They have been utilized successfully in 
modeling various complex spatio-temporal behavior of non-equilibrium dissipative systems such as fluid turbulence, 
anomalous diffusion, and financial markets 

The Levy stable distribution is self-similar to its convolutions. It has a long power-law tail that decays much more 
slowly compared with an exponential decay, which gives rise to infinite variance. However, in practical situations, it 
is usually truncated due to nonlinearity or finiteness of the system. In order to incorporate this fact, Mantegna and 
Stanley || introduced the notion of truncated Levy flight. It is based on a truncated Levy stable distribution with 
a sharp cut-off in its power-law tail. Therefore, the distribution is not self-similar when convoluted, and has finite 
variance. Thus, the truncated Levy flight converges to a Gaussian process due to the central limit theorem in contrast 
to the ordinary Levy stable process. However, as they pointed out, its convergence to a Gaussian is very slow and 
the process exhibits anomalous behavior in a wide range before the convergence. Koponen Q reproduced their result 
analytically using a different type of truncated Levy distribution with a smoother cut-off. 

Dubrulle and Laval || applied their idea to the velocity field of 2D fluid turbulence, and claimed that the truncation 
is essential for the multi-scaling property of the velocity field to appear. They showed that the broken self-similarity 
of the distribution makes a qualitative difference on the scaling property of the corresponding random process, i.e., 
the ordinary Levy stable process exhibits mere single-scaling, while the truncated Levy process exhibits multi-scaling, 
although their analysis was mostly based on numerical calculation and the obtained scaling exponents were rather 
inaccurate. The idea of truncated Levy flights has also been applied to the analysis of financial data such as stock 
market prices and foreign exchange rates MM . Multi-scaling analyses of the financial data have also been attempted 

In this paper, we treat the truncated Levy flights analytically based on the smooth truncation introduced by 
Koponen, and clarify their multi-scaling properties. They exhibit the simplest form of multi-scaling, i.e., bi-fractality, 
due to the characteristic shape of the truncated Levy distribution. Our results may have some relevance to the 
multi-scaling properties of the velocity field in 2D fluid turbulence and of the fluctuation of the stock market prices. 

II. TRUNCATED LEVY DISTRIBUTION 

Let P{x) an probability distribution and its characteristic function, i.e., 

i />oo />oo 

p{x) = — e^e-^dC, e^ (c) = / P{x)e lC - x dx. (1) 

As explained in Feller [|J , (the argument of) the characteristic function for a Levy stable distribution is given by that 
of a compound Poisson process: 

/oo 
(e«« - 1 - Kt{x)) f(x) dx, (2) 
- OO 

where f(x) is a probability distribution of increments and t(x) is a certain centering function. f(x) is assumed to be 

Il{x) ~ \ Cpx- 1 -* (x>0), [6) 
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FIG. 1. Truncated Levy distribution Ptl(x;X) and its n-times convolutions for the symmetric case. The parameters are 
A = 0.001, q = 0.75, and q~p = 0. Corresponding ordinary Levy stable distribution Pl(x) is also shown for comparison. Each 
convoluted distribution is rescaled as n 




FIG. 2. Truncated Levy distribution Ptl(x;X) and its n-times convolutions for the one-sided case. The parameters are 
A = 0.001, a = 0.75, and q — p = —1. Corresponding ordinary Levy stable distribution Pl(x) is also shown for comparison. 
Each convoluted distribution is rescaled as ri 



2 



where C > is a scale constant, 0<a<2,p>0,f/>0, and p + q = 1. The function t(x) is chosen as for 
< a < 1 and x for 1 < a < 2. By integration, we obtain 

MQ=CT(-a)\C\ a (cos^± i ( (Z -p)sin^) (a ^ 1), (4) 

where the upper sign applies when C > 0, and the lower for £ < 0. This is a well-known form of the Levy stable 
characteristic function, and we obtain a Levy stable distribution Pl(x) through an inverse Fourier transform (see 
Figs. and §). 

This characteristic function satisfies jiV'l(C) = ^.t/C^ a C); which means that the corresponding probability distri- 
bution is stable to convolution, i.e., 

Pl{x)=n- 1 ' a P 1 L {n- 1 ^x). (5) 

Here P n {x) denotes a n-times convoluted distribution of P 1 (x) = P(x), i.e., P n (x) = (2tt)^ 1 J e^^'e^^dC,. The 
Levy stable distribution Pl(x) is symmetric when q — p = 0, and one-sided when q — p = ±1 and < a < 1.Q It has 
a power-law tail of the form and the absolute moment (x q ) := \x\ q PL(x)dx does not exist for q > a. 

Now, let us truncate this Levy stable distribution following Koponen Q. We introduce a cut-off parameter A > 
and truncate the original fh(x) in Eq. (||) as 



Cq |x|- 1 - Q e- A l ;c l (x < 0), 
Cp x- l - a e~ Xx (x > 0). 



, / s. i of/ x- e 1 1 i^x ui, 



For the case < a < 1, t(x*) can be omitted, and we obtain by integration 

iPTL(C,\)=Cr(-a){q(\ + t (r+p(\- t (r-\ a }, (7) 
or, by expanding the first two terms 

^tl(C; A) = GT(-a) { (A 2 + C 2 ) a/2 (cos aB + i(q - p) sin off) - A Q } , (8) 

where 9 = arctan (CI A) . Apart from the scale constant, this is the characteristic function of truncated Levy distribution 
given by Koponen.n For the case 1 < a < 2, we use a centering function t(x) = x and obtain 

<Atl(C; A) - CT(-a) {f?(A + *C) tt + p(A - *C) a - A a - iaA^^g - p)C} . (9) 

In this case, an extra term appears in addition to Eq. (0), which induces a drift of the probability distribution when 
(j-p^O. It can easily be seen that in the limit A — > +0, these characteristic functions go back to the Levy stable 
characteristic function Eq. (||). 

We obtain a truncated Levy distribution Ptl(x; A) through an inverse Fourier transform from Eq. (0) or Eq. (|^). 
The parameter A modifies the behavior of i/'tlCC; A) when £ is comparable to A and removes the singularity of the 
characteristic function at the origin. It thus changes the behavior of Ptl(x;X) when x ~ A -1 and introduces an 
exponential cut-off to the power-law tail of Ptl(x; A). Therefore, all absolute moments (x q ) of Ptl(x; A) are finite in 
contrast to the case of the Levy stable distribution Pl(x). 

The characteristic function ^tl(C)A) is infinitely divisible but no longer stable. The convolutions of its corre- 
sponding probability distribution cannot be collapsed by scaling and shifting the variable x. However, they can 
still be collapsed by scaling both x and A appropriately; the characteristic function ^tl(CjA) satisfies mpTL(C,]^) = 
^TL(n 1/a C;n^ a X), which means that the corresponding n-times convoluted probability distribution satisfies 

P^ L (x;X) = n- 1 / a P^ L (n- 1/a x;n 1/a X). (10) 
We utilize this fact for calculating the scaling exponents of the truncated Levy flights. R 



1 The distribution is one-sided to the right (a; > 0) when q — p — —1 and to the left (x < 0) when q — p = 1. 
2 Note the misprint of Eq. (3) in Ref. [§. It should read lnP(fe) = c{c - (k 2 + v 2 y /2 / cos(tt7//2)... } . 

3 In the context of random processes we discuss below, this operation can be considered as a renormalization-group transfor- 
mation. The limiting Levy stable process corresponds to a fixed point, and the exponent 1/a can be regarded as a sort of 
critical exponent. 
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In Figs. [T] and g, truncated Levy distributions Ptl{x; A) and their convolutions are displayed for symmetric (q — p = 
0) and one-sided (q — p = — 1) cases in comparison with the corresponding ordinary Levy stable distributions Pl{x) in 
a log-log scale. As expected, each truncated Levy distribution has a cut-off at x ~ A -1 = I0 3 after a power-law decay. 
The cut-off position gradually approaches the origin with the convolution, and the self-similarity of the convoluted 
distribution is broken in the tail part. 

We can extend the power-law decaying part arbitrarily longer by making A smaller. Hereafter, we assume the cut-off 
parameter A to be very small, i.e., the cut-off is far away from the origin, since we are interested in the transient 
anomalous behavior of the corresponding random process before it converges to a Gaussian due to the central limit 
theorem. 



III. TRUNCATED LEVY FLIGHTS 

The truncated Levy flight (3|-[| is a temporally discrete stochastic process characterized by the truncated Levy 
distribution Ptl(x;X). At each time step i, a particle jumps a random distance x(i) chosen independently from 
Ptl(x; A). The position y(i) of the particle started from the origin is given by y(i) = J2]=i x ti)- Here we consider 
two representative cases of truncated Levy flights, i.e., the symmetric case (q — p = 0, < a < 2) and the one-sided 
case (q — p = —1, < a < 1). 

Figure |3j shows typical realizations of the jump x(i) and the position of the particle y(i) for the symmetric case. 
The time sequence of the jump x(i) is intermittent; x(i) mostly takes small values but sometimes takes very large 
values. Correspondingly, the movement of the particle y{i) is also intermittent. This intermittency gives rise to the 
anomalous scaling property of the trace of y(i) in which we are interested. 

Figure |i| displays typical realizations of x(i) and y(i) for the one-sided case. Since the probability distribution 
Ptl{x] A) vanishes for x < 0, each jump takes a positive value and the position of the particle increases monotonically. 
As in the symmetric case, their time sequences are intermittent. 
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FIG. 3. Typical realization of random sequence x(i) (top) and corresponding truncated Levy flights y(i) (bottom) for 
A = 0.001, a = 0.75, and q — p — (symmetric case). 





IV. MULTI-SCALING PROPERTIES 

In order to characterize the intermittent time sequences shown in Figs. |^ and ^, we introduce multi-scaling analysis. 
It has been employed successfully in characterizing velocity and energy dissipation fields of fluid turbulence and rough 
interfaces of surface growth phenomena Jl(],[ll| . 

Multi-scaling analysis concerns a partition function of the measure defined suitably for the field under consideration. 
Here we focus on the apparent similarity of the intermittent time sequences of truncated Levy flights to those of fluid 
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FIG. 4. Typical realization of random sequence x(i) 
A = 0.001, a — 0.75, and q — p — —1 (one-sided case). 
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turbulence, i.e., the similarity of y(i) in Fig. || to the velocity field, and that of x(i) in Fig. ^ to the energy dissipation 
field. 

Thus, we apply the "multi-affine" analysis to the trace of y(i) for the symmetric sequence. The measure h(n) is 
defined as the absolute hight difference of y(i) between two points separated by a distance n, i.e., h(n) = \y(i+n)—y(i)\. 
The distribution of h(n) does not depend on i, since the increment of this process is statistically stationary. The 

partition function is then defined as Zh(n;q) — (h{n) q ) = ( Y?jJi+i x ti) )> where (...) denotes a statistical average. 
This function is called "structure function" in the context of fluid turbulence. 

On the other hand, for the one-sided case, we focus on the trace of x(i) and apply the "multi-fractal" analysis. 
The measure m(n) is defined as the area below the trace of x(i) between two points separated by a distance n, i.e., 

TO ( n ) — Sj=7+i an d the partition function is defined as Z m (n; q) — (m{n) q ) — x U)j )• 

These partition functions are expected to scale with n as Z^n; q) ~ rfi^ and Z rn (n; q) ~ n T ^ for small n. Further, 
if these scaling exponents £(g) and r(q) exhibit nonlinear dependence on q, the corresponding measures h(n) and m(n) 
are called multi-affine and multi-fractal, respectively, n 

In Fig. ||, we display the partition functions Z m (n; q) for several values of q for the one-sided case in a log-log scale. 
As can be seen from the figure, each partition function exhibits power-law dependence on n for small n. We obtain a 
similar figure for the partition functions Zh(n;q) for the symmetric case. The corresponding scaling exponents £(<?) 
and r(q) are shown in Fig. |^. Each curve exhibits strong non-linear dependence on q; it is linear for small q and 
constant for large q. Thus, the sequences generated by truncated Levy flights clearly possess multi-scaling properties. 



V. DERIVATION OF THE SCALING EXPONENTS 



Now let us derive the scaling exponents ((q) and r(q) from the characteristic function Eq. (0). It is clear from the 
definition of the partition functions that they can be calculated once we know the probability distribution of the sum 
of random variables z(n) := Y]j—i Since the jumps are independent from each other, the probability distribution 
of z(n) is given by a n-times convolution of the truncated Levy distribution Ptl(x; A), i.e., by P^ L (x; A). 



4 The multi- fractal partition function Z m (n;q) defined here is different from the traditional one that is defined as 
N (n) -1 {m(n) q } , where N(n) is a number of boxes of size n that are needed to cover the whole sequence. This makes a 
difference of —I to the scaling exponent r(g) for the one-dimensional case considered here. 
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FIG. 5. Partition functions Z m (n;q) for q 
top one to q — 3.0. 



0.0 - 3.0 at intervals of 0.2. The bottom line corresponds to q — 0.0, and the 



As we explained previously, Pf L (x; A) can easily be obtained from the original Ptl(x; A) by scaling x and A. Making 
use of this fact, the <?-th absolute moment (z(n) q ) of z(n) can be calculated as 



{z(n) q ) x =A z q P£ L (z;\)dz = A z q n~ 1 l a P TL (rr 1 ^ 'z;n^ a \)dz 



*' a A / z q P TL (z; n 1 / a X)dz = n q ' a (z(l) q ) nl/aX , 



(11) 



where the constant A is 2 for the symmetric case and 1 for the one-sided case. Here we explicitly indicated the 
parameter A of the distribution Ptl(x; A) as the subscript of the average. Note that if {z(l) 9 )\ does not depend on 
A, the scaling exponent is simply given by a linear function q/a, and the process does not exhibit multi-scaling. This 
is the case for the ordinary Levy stable distribution. 

Thus, all we need to calculate is the moment (z(l) q )\. However, of course, an analytical expression for Ptl(x; A) is 
not attainable except a few specific cases. Here we adopt an approximation which utilizes the facts that the truncated 
Levy distribution is different from the ordinary Levy stable distribution only in the tail part, and that it has a power- 
law decaying part with an exponent —1 — a in the middle (see Figs. |l| and |^). Therefore, it is expected that the 
moment whose degree q is lower than a is almost the same as that obtained from the ordinary Levy stable distribution, 
and the moment for q > a is mostly determined by the asymptotic tail of the truncated Levy distribution. Note that 
the moment for q > a does not exist for the ordinary Levy stable distribution. 

(i) Lower moments (0 < q < a). 

We can approximate Ptl{x; A) by Pl{x) in this case, since they are different only in their tail parts, which does 
not contribute to the moments lower than a significantly. Therefore, (z(l) q )\ does not depend on A approximately, 
and we obtain (z(n) q ) ~ const. n q l a . Thus, the scaling exponents and r(q) are given by q/a for < q < a. The 
broken self-similarity of the distribution is not important in this regime. 

(ii) Higher moments (q > a). 

Since the tail of the distribution mainly contributes to the higher moments, we can calculate (z(l) q )\ approximately 
if we know the asymptotic form of the distribution. For this purpose, we expand the characteristic function as in the 
case of the series expansion of ordinary Levy stable distribution [El. By expanding the integrand, the truncated Levy 
distribution can be expressed as 



P TL (x;X) = e- cr (-^ a l- 
Zn 



exp [CT(-a) {q(X + i() a + p(A - i() a }} e^dQ 
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FIG. 6. Mult i- scaling exponents £(q) and r(q) for the symmetric (top) and one-sided (bottom) cases. The cut-off parameter 
A is set at 0.00001. 
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for the case < a < 1. At the lowest order (k = 1), we recover the original distribution of the increments: 



P TL {x; A) ~ e- cr '- a ' r Cpi- 1 - a e- Al {x > 1). (13) 
With this approximation, the moment (z(l) q )\ is calculated as 

/>oo 

(z(l) 9 } A ~ e- cv{ ~ a)x " Cp / z^^e-^dz = e- cr( - a)xa Cpr(l + q)\- q+a , (14) 



and we obtain 

(z{n) q ) x = ni/ a {z{iy) nl/ax ~ n«/ Q (n^n) ?+Q - n 1 . (15) 

Thus, the scaling exponents Q{q) and r(q) are given by 1 for q > a. 

In summary, we approximately derived the following expression for the multi-scaling exponents C(q) and T(q): 

{q/a (0 < q < a), 
(16) 
1 (q>a). 

Note that the above approximation becomes more and more accurate as we extend the power-law decaying part by 
decreasing A, and this result is exact in the asymptotic limit. In Fig. ^, these theoretical curves are compared with 
the experimental results. Except for small deviations near the transition points, the theoretical results well reproduce 
the experimental results. Although our estimation here is done for the case < a < 1, the theoretical result Eq. jl^ ) 
seems to be also applicable for 1 < a < 2.Q This type of simple multi-scaling is sometimes called "bi-fractality" and 
is known, for example, in randomly forced Burgers' equation |10|. 



VI. DISCUSSION 



We analyzed the multi-scaling properties of the truncated Levy flights based on the smooth truncation introduced 
by Koponen. As Dubrulle and Laval claimed, truncation is essential for the multi-scaling properties to appear. We 
clarified this fact and derived the functional form of the scaling exponents for both symmetric and one-sided cases. 

As we mentioned previously, the cut-off parameter A may represent the finiteness of the system under consideration. 
Then it would be natural to assume A as a decreasing function of the system size L, for example A = L , and we can 
think of the system-size dependence of the truncated Levy flight. Of course, distribution functions Ptl for different 
system sizes cannot be collapsed simply by rescaling x, i.e., finite-size scaling does not hold.[] However, approximate 
finite-size scaling relations still hold separately, if we divide Ptl into two parts at the power-law decaying part, i.e., 
into the self-similar part and the tail part. As we explained, the self-similar part is insensitive to A = L , and Ptl 
for different system sizes collapse to a single curve in this region without rescaling. On the other hand, the tail part is 
asymptotically given by P TL (x;£ _1 ) ~ x"~ 1 ~ a e- x / L , which can be rescaled as L 1+€ *Ptl(xL; L x ) to give a universal 
curve. Thus, we have two different approximate finite-size scaling relations in different regions of x, which is separated 
by the power-law decaying part. Of course, this is closely related to the bi-fractal behavior of the scaling exponents. 
Similar asymptotic finite-size scaling is also reported in the sandpile models of self-organized criticality p3|. 

Recently, Chechkin and Gonchar |l4| discussed the finite sample number effect on the scaling properties of a stable 
Levy process. (They treated only the symmetric case using a different argument from ours, which was more qualitative 
and therefore more general in a sense.) They claimed that, due to the finite sample number effect, "spurious" multi- 
affinity is observed in the numerical simulation and derived the "spurious" multi-scaling exponent. Interestingly, or 
in some sense obviously, their "spurious" multi-scaling exponent is the same as our multi-scaling exponent, since the 
truncation of the power-law by A can also be interpreted as mimicking the finite sample number effect of experiments. 



5 This implies that Eq. (^2|) still has its meaning as an asymptotic expansion for 1 < a < 2. This is proved for the case of 
ordinary Levy stable distribution Jli|. 

6 We may also take a different viewpoint, where A is another tunable parameter independent of the system size L. Then Eq. ( |ic| ) 
may be viewed as a finite-size scaling relation between systems of size L/n and L, similar to that in statistical mechanics. 
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Our calculation presented in this paper is similar to our previous work |15[| on the multi-scaling properties of the 
amplitude and difference fields of anomalous spatio-temporal chaos found in systems of non-locally coupled oscillators. 
The distribution treated there was not the truncated Levy type but had a form like (1 + x 2 ) a / 2 e~ x ^ . Since this 
form is easily generated by a simple multiplicative stochastic process, some attempts have been made to model the 
economic activity using this type of distribution [|l6],[l7| . 

The bi-fractal behavior of the scaling exponent is the simplest case of multi-scaling, while experimentally observed 
scaling exponents usually exhibit more complex behavior. Actually, it has long been discussed in the context of fluid 
turbulence what the shape of the distribution should be to reproduce the experimentally observed scaling exponent 
p0| , pl]p8| . In order to reproduce the behavior of the scaling exponent more realistically in the framework of the 
truncated Levy flights discussed here, introduction of correlations to the random variables will be necessary. Studies 
in this direction are expected. 
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